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Description 
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Nomenclature  (Continued) 


temperature 
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mass  fraction  of  element  i 
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c 

P 
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s 

w 

o 

oo 
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Description 

combustion  products 
particle  surface 
reference  state 
solid 

room  surfaces 

initial  condition 

ambient  state  for  concentrations 

ambient  state  for  temperature 


Summary 


The  results  of  a  brief  preliminary  study  of  the  combustion 
properties  of  carbon  slurry  fuels  are  reported.  The  carbon  slurry 
fuels  were  formulated  by  mixing  noncatalyzed  and/or  catalyzed  medium 
thermal  carbon  black  (roughly  50  percent  by  weight)  with  JP-10. 

The  combustion  of  individual  supported  slurry  drops,  400-1000  pm 
initial  diameter,  was  observed  in  an  open  turbulent  diffusion  flame 
at  atmospheric  pressure.  The  environment  of  the  drops  within  the 
flame  was  known  from  measurements  of  mean  velocity,  temperature  and 
composition.  The  drops  were  supported  on  the  junction  of  a  thermo¬ 
couple.  Shadowgraph  motion  pictures  were  also  obtained  for  the 
process.  This  data  provided  the  variation  of  particle  temperature 
and  diameter  as  a  function  of  time. 

A  two-stage  combustion  process  was  observed.  The  first  stage 
involved  heat-up  and  evaporation  of  the  liquid,  leaving  a  solid 
agglomerate  of  the  carbon  particles  in  the  slurry.  The  second 
stage  involved  heat-up  and  reaction  or  quenching  of  the  agglomerate. 

A  model  was  constructed  for  the  second  stage  of  Lho  process.  This 
involved  analysis  of  heat  and  mass  transfer  processes  in  the  convective 
environment  surrounding  the  agglomerate,  and  the  chemical  reaction  oi 
carbon  with  oxygen,  carbon  dioxide,  and  water  vapor  at  the  surface 
of  the  particle.  Radiative  heat  transfer  between  the  agglomerate  and 
its  surroundings  was  also  considered.  The  predictions  of  the  model 
were  compared  with  measurements  in  the  lean  portions  of  the  flame 
(equivalence  ratios  in  the  range  0.272-0.778). 

The  major  findings  of  the  study  can  be  summarized  as  follows: 

1.  For  present  test  conditions,  all  the  solid  carbon  in  the 
slurry  remained  in  the  agglomerate  and  heat-up  and  reaction 
of  the  agglomerate  dominated  the  lifetime  of  the  particle — 
usually  requiring  90-95  percent  of  tiie  particle  lifetime, 
even  in  locations  where  agglomerate  reaction  was  nearly 
diffusion  controlled. 

2.  While  liquid  heat-up  and  gasification  occurs  in  all  regions 
of  the  flame,  agglomerate  reaction  was  limited  to  regions 
with  equivalence  ratios  in  the  range  0.37-4.0  for  noncatalyzed 
slurries  and  0.21-4.0  for  catalyzed  slurries.  The  carbon 
reaction  was  quenched  at  other  locations. 

3.  The  model  yielded  good  predictions  of  agglomerate  reaction 
rates  in  the  range  tested  (equivalence  ratios  of  0.272-0.778, 
where  the  main  reactant  ;s  oxvgen).  Conditions  near  both  the 
diffusion  and  reaction  rate  controlled  limits  were  included 

in  the  comparison.  This  agreement  was  obtained  by  selecting  a 
single  empirical,  area/reactivity  multiplication  factor,  which 
allows  for  the  pore  and  surface  structure  of  the  carbon  ag¬ 
glomerate  and  the  effect  of  a  catalyst. 


4.  The  model  yielded  good  predictions  of  agglomerate-life- 
histories  over  the  same  range  of  conditions.  Reasonable 
predictions  of  heat-up  times  and  particle  temperatures 
indicate  that  the  thermal  aspects  of  the  model  are  adequate 
over  the  present  test  range. 

5.  Adding  catalyst  increased  agglomerate  burning  rates,  for 
conditions  not  too  near  the  diffusion  controlled  limit. 

This  effect  could  be  treated  in  the  theory  by  selecting 
a  higher  area/reactivity  multiplication  factor  for  the 
catalyzed  slurry.  The  presence  of  catalyst  also  extended 
the  lean  limit  for  quenching. 

6.  Some  particle  shattering  was  observed  with  catalyzed 
slurries  in  the  lean  region  of  the  flame. 

These  conclusions  are  based  on  test  results  with  relatively  large 
particles,  400-1000  Um  in  diameter.  In  particular,  the  effect  of  pore 
diffusion  is  likely  to  vary  with  particle  size.  Therefore,  the  area / 
reactivity  multiplication  factors  selected  to  correlate  the  present 
data  are  not  considered  to  be  appropriate  for  other  test  conditions, 
without  further  testing.  The  effect  of  pressure  and  combustion 
properties  in  the  fuel  rich  region  (equivalence  ratios  of  0.778-4.0) 
were  also  not  assessed  in  this  study.  This  implies  considerable 
uncertainty  concerning  agglomerate  reaction  with  carbon  dioxide, 
water  vapor,  and  other  species  present  at  high  equivalence  ratios.  In 
order  to  remove  these  uncertainties,  it  is  recommended  that  additional 
work  be  undertaken  using  particle  sizes  more  representative  of  com¬ 
bustion  chamber  conditions,  as  well  as  test  environments  extending 
to  richer  mixture  ratios. 


1.  Introduction 


1 . 1  General  Objectives 

The  propulsion  system  of  airbreathing  missiles  is  frequently 
volume  limited.  This  implies  that  the  volume  of  the  propulsion  system 
is  a  critical  factor  in  determining  missile  range  and  performance — as 
opposed  to  propulsion  system  weight.  Since  fuel  volume  comprises  a 
significant  fraction  of  propulsion  system  volume,  fuels  that  have  a 
large  chemical  energy  release  per  unit  volume  (high  volumetric  energy 
density)  are  particularly  attractive  for  such  applications. 

There  has  been  continuous  development  of  hign  volumetric  energy 
density  liquid  fuels  to  meet  these  needs  [1].  However,  it  has  become 
increasingly  difficult  to  achieve  further  improvements  using  liquids 
alone.  This  has  resulted  in  increased  consideration  of  slurries, 
since  these  fuels  can  be  handled  and  burned  similar  to  liquids  and 
yet  have  increased  density  due  to  the  presence  of  solids  [2-5].  The 
objective  of  the  present  investigation  was  to  conduct  a  preliminary 
investigation  of  the  combustion  properties  of  carbon  slurry  fuels  of 
interest  to  the  Air  Force  [2,3]. 

The  carbon  slurry  fuels  considered  in  this  investigation  were  tormu- 
lated  by  Suntech,  Inc.,  Marcus  Hook,  PA,  for  use  during  combustor  tests 
by  the  AiResearch  Manufacturing  Company  of  Arizona  [2,3].  While  these 
fuels  are  representative  of  candidate  carbon  slurry  fuels,  they  do  not 
necessarily  represent  the  most  attractive  formulations  since  fuel  de¬ 
velopment  efforts  are  still  in  progress.  Nevertheless,  tiie  materials 
should  still  provide  a  reasonable  indication  of  carbon  slurry  combustion 
properties,  useful  for  subsequent  fuel  development  efforts. 

The  properties  of  the  fuel  considered  during  the  present  investi¬ 
gation  are  summarized  in  Table  1.  The  fuels  consist  of  a  medium  tuermal 
carbon  black,  with  JP-10  as  the  liquid  fuel.  Two  slurries  were  tested, 
each  containing  roughly  50%  carbon  by  weight.  The  ultimate  carbon 
particle  size  was  0.3  pm.  One  slurry  was  catalyzed  with  a  proprietary 
lead  compound . 

The  overall  combustion  characteristics  of  carbon  slurry  fuels  nave 
been  examined  in  gas  turbine  combustors  [3],  and  in  well-stirred 
reactors  [5].  The  objective  of  the  present  investigation  was  to  examine 
slurry  combustion  on  a  more  fundamental  level,  by  observing  the  com¬ 
bustion  of  individual  slurry  drops.  The  experimental  technique  involved 
observations  of  single  slurry  drops  (400-1000  pm  in  diameter).  The 
drops  were  supported  at  various  positions  within  a  turbulent  diffusion 
flame  in  order  to  provide  an  environment  representative  of  actual  com¬ 
bustion  chamber  conditions.  Measurments  were  made  of  the  variation  ol 
particle  diameter  and  temperature  as  a  function  of  time.  The  process  was 
also  modeled  theoretically  and  predictions  and  measurements  were  compared. 
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Table  1 


Properties  of  the  Carbon  Slurry  Fuels 


a 


Noncatalyzed 

Catalyzed13 

Designation 

790-928 

790-942 

Liquid 

JP-10 

JP-10 

Q 

Dispersed  Carbon  (wt  %) 

50.4 

49.2 

aFuels  supplied  by  R.  S.  Stearns,  Suntech,  Inc., 

P.  0.  Box  1135,  Marcus  Hook,  PA  19061. 

The  catalyst  is  a  proprietary  lead  compound  [6], 

CThe  dispersed  carbon  is  a  medium  thermal  carbon  (carbon 
black)  having  an  ultimate  particle  size  of  0.3  pm. 
Further  properties  of  this  material  are  provided  in 
Ref.  2. 
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1.2  Previous  Studies 


The  present  study  is  an  extension  oi  earlier  work  in  this 
laboratory  [3,4].  The  test  fuel  and  many  aspects  ol  the  experiments 
are  similar.  The  general  test  procedure  involved  observations  of 
slurry  drops  supported  in  an  open  turbulent  diffusion  flame. 

The  major  findings  of  Ref.  4  are  as  follows: 

1.  The  combustion  of  a  slurry  drop  occurred  in  four  stages: 

1)  heat-up  of  tne  drop  and,  at  times,  ignition  of  liquid 
fuel  vapors;  2)  evaporation  or  combustion  of  the  liquid  fuel 
accompanied  by  agglomeration  of  the  carbon  particles  in  the 
slurry;  3)  further  heating  of  the  carbon  agglomerate  once  all 
the  liquid  has  evaporated;  and  4)  fragmentation,  reaction  or 
quenching  of  the  agglomerate. 

2.  Combustion  of  the  carbon  agglomerate  only  occurred  over  a 
limited  equivalence  ratio  range  (0.21-4  for  the  ratal  wed 
slurry,  0.37-4  for  the  noncatalyzed  slurry). 

3.  Even  in  regions  where  the  carbon  agglomerate  reacted,  Lite 
heat-up  and  reaction  of  Lite  agglomerate  was  the  slowest  step 
in  the  combustion  process,  generally  requiring  20  times 
longer  than  the  time  required  to  gasify  the  liquid. 

Figure  1  is  a  scanning  electron  microscope  photograph  of  t  he 
surface  of  a  partly  reacted  agglomerate.  The  light  areas  on  Lite 
photograph  are  the  surface  of  the  solid  agglomerate,  the  dark  areas 
represent  pores.  The  agglomerate  is  clearly  very  porous,  with  the 
carbon  surface  area  available  for  reaction  much  greater  than  the 
apparent  outer  surface  area  of  the  agglomerate. 

These  general  observations  concerning  the  combustion  properties 
of  carbon  agglomerates  are  similar  to  results  found  for  coal  slurries 
[7,8].  In  the  case  of  coal,  however,  the  presence  of  ash  and  volatil 
material  in  the  coal,  materially  influences  the  details  of  the  proves 

1 . 3  Specific  Objectives 

Previous  investigations  have  demonstrated  the  i'-'port  ain-i  ot 
agglomerate  formation  and  reaction  duri.ig  the  combustion  ot  carbon 
slurries.  Some  information  concerning  the  regions  of  t  lames  where 
agglomerates  react  lias  also  been  generated.  however,  exist  in  data  o 
agglomerate  reaction  rates  are  very  limited.  Furthermore,  no  attempt 
to  model  the  agglomerate  reaction  process  lias  been  reported. 

With  this  status  in  mind,  this  investigation  extended  the  result 
of  Refs.  3  and  4.  employing  a  similar  experimental  arrangement.  I u< 
specific  objectives  of  the  study  were  as  lollows: 
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i.  Complete  additional  measurements  ot  the  variation  oi 
particle  size  and  temperature  witli  time  lor  drops  sup¬ 
ported  at  various  positions  within  a  turbulent  diliusion 
flame.  In  particular,  generate  more  complete  information 
on  agglomerate  burning  rates  as  a  1  unction  of  position  in 
tiie  flame. 

ii.  formulate  a  model  tor  the  reaction  of  the  carbon  agglom¬ 
erate,  using  existing  inlormation  on  carbon  reaction 
rates.  Fmploy  the  model  to  compute  agglomerate-1 i fe- 
histories  and  compare  predictions  and  measurements. 

The  supported  particles  used  during  Lliis  investigation  had  diameters 
in  the  range  400-1000  pm,  which  is  much  larger  than  typical  combustor 
applications.  For  preliminary  work,  however,  Lhis  approach  had  ad¬ 
vantages  since  more  detailed  measurements  can  be  made  with  large 
particles,  helping  to  provide  more  insight  concerning  shirr'  comnustiou 
properties . 

2.  Test  Arrangement 

1 . 1  Apparatus 

The  test  apparatus  involved  modification  ol  an  earlier  arranvemeiit 
employed  for  a  studs  of  gas  and  spins  diffusion  flames  |‘f,lu|.  A  skclNi 
ot  the  modified  apparatus  appears  in  Fig.  .’ .  i  he  arrangement  consist., 
an  injector  flowing  propane  gas,  winch  is  burned  as  a  diffusion  flame 
in  stagnant  room  air.  The  test  drop  is  supported  from  either  a  tnermo- 
couple  or  a  quartz  probe  at  various  locations  in  the  flame.  1  he  drops 
are  also  photographed  with  a  motion  picture  camera,  as  shadowgraphs . 

The  propane  gas  injectoi  wa->  posit  ioned  mat  the  hot  tom  ol  the  lest 
stand  and  oriented  vertical  Is  upward.  Lhc  ten  stand  was  an  aie.i  1..’  mm 
square  and  i  m  high  enclosed  wiln  a  single  thickness  ol  lh  mesh  s<  reen. 

1'he  test  stand  was  located  within  a  room  Having  d  i  mens  i  oils  4m  x  7m  x  4m 
high.  Combustion  products  were  removed  turough  an  exhaust  fan  located 
near  the  ceiling  ot  the  test  cell.  1  he  injector  was  seated  on  a  support 
which  allowed  a  vertical  movement  ol  1.2  m.  Iwo  traversing  mecii.m  i  sms 
were  used  to  position  the  nozzle  in  the  Horizontal  plane. 

The  nozzle  was  a  lull-cone  air  atomizing  in  je.  lor,  with  no  swirl, 
manufactured  by  the  Spraying  Systems  i  ompans  (Model  I/s  1  i  luid 

nozzle  and  Model  7 1 4  7  air  nozzle).  Ihe  outlet  diameter  m  tm  nozzle 
was  1.144  mm. 

The  gas  flame  was  stabilized  near  the  exit  of  the  iujeiloi  using  an 
array  of  four  hydrogen  capillarv  flames.  flie  Name  tubes  wen  mounted 
in  a  symmetrically  opposed  pattern  injecting  toward  the  centerline  ol 
the  injector  passage.  The  flow  passages  tor  the  hydrogen  were  slots 
(J .  4  mm  wide  and  2  mm  long,  witli  the  long  axis  parallel  to  the  injector 
centerline.  The  lower  edges  ot  the  slots  were  1 . H  mm  above  the  lace  of 
the  injector.  ihe  exit  planes  ol  the  slots  were  h . 4  mm  apart . 


The  propane  gas  flow  rate  was  metered  with  a  Matheson  Model  604 
rotameter  and  controlled  by  a  Harris  Model  2515  pressure  regulator 
with  an  output  capacity  of  0-0.69  MPa.  The  flow  rate  of  the  Hydrogen 
gas  was  metered  with  a  Matheson  Model  601  rotometer  and  controlled 
with  a  Matheson  Model  111  pressure  regulator  having  an  0-1.38  MPa 
output  capacity.  The  rotameters  for  both  the  propane  and  hydrogen 
flows  were  calibrated  with  a  Precision  Scientific  Company  wet-test 
meter  (2.83  ml/rev). 

The  fuel  drops  were  mounted  with  a  hypodermic  needle.  The  flame 
was  deflected  from  the  mount  region  until  the  droplet  was  in  place. 

The  deflector  was  constructed  from  a  200  x  305  mm  sheet  of  stainless 
steel,  having  a  thickness  of  3  mm.  The  drop  combustion  process  was 
initiated  by  removing  the  deflector  with  a  pneumatic  cylinder  arrange¬ 
ment. 

2.2  instrumentation 


2.2.1  Drop  Environment  Measurements 

The  drop  environment  at  each  test  location  within  the 
flame  was  measured.  This  involved  measurements  of  mean  velocity, 
temperature  and  composition.  Mean  velocities  were  measured  with  a 
laser-Doppler  anemometer,  which  employs  frequency  shifting  so  that 
flow  reversals  could  be  detected.  Mean  temperature  was  measured  with 
a  50  pm  diameter  platinum/platinum-10%  rhodium  thermocouple  (with  a 
radiation  error  of  less  than  35  C  over  the  test  range).  Mean  concen¬ 
trations  were  measured  by  extracting  gas  samples  isokinetically  and 
analyzing  them  with  a  gas  chromatograph .  The  details  of  these 
techniques  are  presented  elsewhere  [9,10]. 

2.2.2  Drop  Measurements 

The  drop  combustion  process  was  observed  with  a  16  nun, 
Photosonics,  Model  16-B,  motion  picture  camera.  The  camera  optics 
provided  a  2:1  magnification.  Tiie  camera  was  powered  with  a  Kepco, 

SM  36-5  AM  d.c.  power  supply.  The  film  speed  was  indicated  with  a 
timing  light  on  the  camera,  activated  with  an  Adtrol  Electronic  pulse 
generator,  Model  501.  Kodak  Plus-X  reversal  1 ilm  was  used  tor  all 
tests. 

Backlighting  for  shadowgraph  measurements  was  provided  by  a  Pek, 
Model  401A,  75W  mercury  arc  lamp.  The  light  Iron)  the  arc  was  collimated 
and  directed  toward  the  drop  location.  A  dil laser  screen  was  employed 
behind  the  drop  to  equalize  the  1 ight  intensity  ot  the  background.  I  he 
background  intensity  was  adjusted  so  that  envelope  I  lames  around  the 
drop  could  also  be  observed. 

The  film  records  were  analyzed  on  a  frame  by  Irame  basis,  using  a 
Vanguard  motion  picture  analyzer.  Photographs  ol  objects  ol  known  size 
at  the  drop  location  provided  a  calibration  of  distances  on  the  t  ilm. 
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Particle  temperatures  were  measured  by  mounting  the  test  drop  at 
the  junction  of  a  thermocouple.  A  Pt/Pt-10%  Rh  thermocouple  was 
used,  constructed  of  75  pm  diameter  wires.  In  order  to  help  support 
the  drop,  the  junction  was  placed  within  a  round  bead,  formed  with 
Sauereisen  cement  (300-400  pm  in  diameter).  The  output  of  the 
particle  thermocouple  was  recorded  digitally  on  a  Nicolet  Explorer 
Oscilloscope,  Model  206,  and  reduced  to  temperature-time  data  on  a 
computer . 

2.3  Test  Conditions 


Tests  were  conducted  using  the  burner  operating  condition 
summarized  in  Table  2.  The  drops  were  mounted  along  the  centerline  of 
the  flow.  Relatively  complete  structure  measurements  have  been  reported 
for  this  flame,  as  noted  earlier  [9].  Details  of  interest  to  this  in¬ 
vestigation  will  be  discussed  later. 


3.  Theoretical  Considerations 


3 . 1  Description  of  the  Model 

Results  obtained  during  this  and  earlier  investigations  indi¬ 
cate  that  slurry  drop  combustion  is  primarily  a  two-stage  process  in¬ 
volving  gasification  of  liquid  in  the  first  stage  and  heat-up  and 
reaction  of  the  agglomerate  in  the  second  stage.  Since  the  agglomerate 
does  not  actively  participate  in  the  first  stage,  other  than  influencing 
the  density  of  the  particle,  it  is  convenient  to  treat  these  two  stages 
separately . 

The  following  analysis  will  be  limited  to  the  second  stage  —  involving 
heat-up  and  reaction  of  the  agglomerate.  The  present  investigation 
emphasized  this  aspect  of  the  problem  since  reaction  of  the  agglomerate 
is  the  slowest  step  in  the  process.  Therefore,  agglomerate  reaction 
generally  controls  combustor  performance.  Furthermore,  agglomerate 
reaction  presents  the  greatest  uncertainties;  in  contrast,  gasification 
of  liquid  drops  is  reasonably  well  understood  [11]. 

The  initial  condition  of  the  present  analysis  is  the  point  where 
all  the  liquid  has  vaporized  and  the  agglomerate  is  at  the  wet-bulb 
temperature  appropriate  for  liquid  vaporization  (the  wet-bulb  tempera¬ 
ture  is  slightly  below  the  boiling  temperature  of  the  liquid  lor  a  pure 
fuel  like  JP-10  [11]).  The  objective  of  the  analysis  is  to  predict  the 
subsequent  heat-up  and  reaction  of  the  agglomerate.  This  involves 
determining  the  variation  of  particle  temperature  and  diameter  as  a 
function  of  time. 

3 . 2  Agglomerate  Reaction  Model 

Since  the  test  droplets  were  supported  at  fixed  locations, 
the  motion  of  the  particle  is  not  considered.  The  lormulation  was 
limited  to  determining  the  rate  of  consumption  of  carbon  in  the 
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Table  2 


Burner 


Flame  Characteristics 


a 


Fuel: 

propane 

Fuel  flow  rate: 

176  mg/s 

Initial  fuel  jet  velocity: 

38.7  m/ s 

Initial  fuel  jet  diameter: 

1 . 194  mm 

Jet  Reynolds  number: 

23600 

Orientation : 

vertical  (upwar. 

Flame  Height: 

460  mm  (visual) 

Injector  Thrust: 

15.6  mN 

Hydrogen  flow  rate: 

0 . 126  mg/ s 

Ambient  and  injector  inlet  temperature: 

296  K 

Ambient  pressure: 

97  kPa 

aNozzle  is  Spraying  Systems  Co.  Model  1/4  J  2050 
Model  67147  air  nozzle,  air  atomizing  injector. 

fluid  nozzle  and 

agglomerate,  the  rate  of  heat  transfer  to  the  particle,  and  the  rate 
of  chemical  energy  release  at  the  surface  of  the  particle. 

A  sketch  of  the  model  appears  in  Fig.  3.  The  major  assumptions 
of  the  analysis  are  as  follows: 

1.  The  particle  is  assumed  to  be  spherical  and  particle  spacing 
is  large  so  that  the  surroundings  can  be  taken  as  an  infinite 
medium  having  uniform  properties. 

2.  The  gas  flow  fieiu  around  the  particle  is  assumed  to  be 
quasisteady,  i.e.,  the  flow  field  has  properties  equivalent 

to  a  steady  state  flow  field,  for  the  same  boundary  conditions, 
at  each  instant  of  time. 

3.  Flow  fields  around  spheres  are  too  complex  to  analyze  exactly; 
therefore,  the  conventional  film  theory  approximation  is 
employed  to  represent  the  effect  of  convection  [11] .  This 
assumption  implies  that  gas  phase  transport  processes  can  be 
represented  as  a  spherically  symmetric  flow  through  a  stagnant 
layer  whose  outer  boundary  is  fixed  empirically. 

4.  The  particle  is  assumed  to  consist  of  pure  carbon  and  lias  a 
uniform  temperature  at  each  instant  of  time. 

5.  Solid  carbon  is  only  assumed  to  react  at  the  apparent  surface 
of  the  particle.  The  effect  of  pores  and  reaction  within  the 
particle  is  treated  by  introducing  an  empirical  area/react  ivi ly 
multiplication  factor.  Reaction  of  solid  carbon  with  oxygen, 
carbon  dioxide, and  water  vapor  is  considered. 

6.  The  gases  surrounding  the  particle  are  assumed  to  be  represented 
by  their  mean  properties  and  the  effect  of  turbulent  fluctuations 
is  ignored.  The  actual  gas  mixture  was  approximated  by  oxygen, 
nitrogen,  carbon  dioxide,  carbon  monoxide,  water  vapor,  hydrogen 
and  hydroxyl  radicals — in  thermodynamic  equilibrium. 

7.  Only  concentration  diffusion  is  considered,  employing  an 
effective  binary  diffusion  coefficient  for  all  diffusing  species. 
All  species  are  assumed  to  have  equal  specific  heats  (assumed 
constant)  and  molecular  weights,  for  the  purpose  of  computing 
diffusion  and  energy  transport  rates  in  the  flow.  Constant 
average  values  are  used  for  all  properties  in  the  gas  phase. 

8.  The  gas  is  assumed  to  be  transparent  to  radiation,  however, 
provision  is  made  for  particle  radiation  to  the  confinement 
surrounding  the  flow. 

9.  The  pressure  within  the  flow  field  is  assumed  to  be  constant. 

10.  The  gases  maintain  thermodynamic  equilibrium  at  the  particle 
surface  however,  the  flow  is  chemically  frozen  elsewhere. 
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These  assumptions  are  generally  similar  to  those  employed  in 
recent  analyses  of  carbon  particle  reaction  [12,13]. 


3.2.1  Gas  Phase  Transport 

The  governing  equations  consist  of  equations  for  con¬ 
servation  of  mass,  energy  and  species.  The  latter  equation  can  be 
represented  more  conveniently  in  the  present  analysis,  where  equal 
binary  dif f usivities  are  assumed,  as  equations  for  conservation  of 
elements — since  elements  are  conserved  during  chemical  reactions. 
Employing  the  assumptions  of  the  analysis,  the  following  equations 
apply  [14]. 


2 

pvr  = 


(1) 


,,  2  dh 

m  r  — 
P  P  dr 


_A_  d_ 

C  dr 
P 


,  2dh, 
(r  dr) 


(2) 


„  2dYi 

n  r  _ — 

p  p  dr 


d  2  d  % 

=  ^ ^Yi»’  1  ^ 


(3) 


Four  element  fractions  are  considered  in  the  analysis,  Y^  ,  Y^  ,  Y^  and'1)1  • 

Seven  species  were  considered:  0  ,  N  ,  C0„,  CO,  H^O,  and  OH.  following 
Refs.  12  and  13,  0  was  ignored  on  the  grounds  that  its  concentration  is  too 
small,  for  temperatures  below  4000K,  for  it  to  contribute  significantly 
to  the  reaction  of  carbon.  Species  and  element  concentrations  are  related 
by  the  following  mass  balances 


Y. 

l 


V.  .  Y. 
ij  3 


i  =  1,4 


(4) 


where  V  is  the  mass  fraction  of  the  ith  element  in  the  jth  species. 

ij 


di 


wMflMatfttilMii 


The  enthalpy  appearing  in  Eq.  (2)  is  given  by 


h  = 


i=l 


Y.h. 
1  1 


(5) 


where 


h. 

i 


(T-Tref)  +hf. 

i 


ref 


The  boundary  conditions  for  Eqs.  (2)  and  (3)  are  as  follows: 


r=r  ,  h=h  ,  ^.=^.  ;  r=r  ,  Y.=Y.  ;  r=r  ,  h=n  (7) 

p  Dll  co  *11  oj  ’  .jo 

P  «>  T 


The  ambient  conditions  are  specified.  The  conditions  at  the  particle 
surface  must:  be  determined  from  the  analysis  of  surface  reaction.  This 
will  be  discussed  later. 

The  general  case  is  treated  in  Eq.  (7),  where  the  film  radii  for 
mass  and  heat  transfer,  r^  and  r^  ,  are  not  the  same.  These  radii  are 

found  employing  the  usual  correlations  from  film  theory  [11].  For  mass 
transfer. 


r  /r  -  Sh/  (S'n-2)  (8) 

CO  p 

while  for  heat  transfer 

r  /r  =  Nu/ (Nu-2)  (9) 

P 


The  Sherwood  and  Nusselt  number  correlations  employed  in  Eqs.  (8)  and 
(9)  are  found  from  the  following  correlation  [11]. 


(Nu  or  Sh)-2 


0.555  Re1/2(Pr  or  Sc)1'3 


;i+1.232/(Re(Pr  or  Sc  )  ^3) )  1^*‘ 


(10) 


Equations  (2)  and  (3)  can  be  integrated,  subject  to  the  boundary 
conditions  of  Eq .  (4),  to  yield 


% 

Y. 


1 


exp(-K/£)  +  Y,p  exp(-K/6)  -  Yioo  exp(-K) 
exp(-K/d)-exp(-K) 


ai; 


(hoo_hp)  exP  (~K/ t,Le)  +  h  exp (-K/yLe)  -  h^exp (-K/Le)  (12) 

h  =  — - *- - - - — - 1 — — _ _ 

exp(-K/(Le)  -exp(-K/Le) 

where 

C  =  r/r  (13) 

6  =  /rp  (14) 

Y  =  roo  /r  (15) 

T  H 


and  K  is  a  nondimensional  particle  mass  loss  parameter,  defined  as 


K  =  m"  r  /pD 
P  P 


(16) 


A  relationship  between  K  and  the  can  be  developed  by  noting 

that  there  is  no  net  mass  flux  of  the  elements  oxygen,  hydrogen  or 
nitrogen  at  the  particle  surface.  Furthermore,  the  net  mass  flux  of 
carbon  at  the  surface  is  equal  to  m",  since  the  particle  is  assumed  to 
be  pure  carbon.  Therefore,  at  £=l,^we  have 


d£ 


i  =  1,4 


i  *  C 


(17) 
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(IB) 


-K(l-Yc) 


Taking  the  derivative  of  Eq.  (11)  and  substituting  into  Eqs.  (17)  and 
(18)  yields  the  following 


Y.  =  Y.  exp (K(6  1-1))  ,  i  -  1,4  ,  i  t  C  (19) 


Yc  =  1-(1-YC  )  exp (K(6  1-1))  (20) 

D  oo 


Employing  the  usual  sign  convention  that  convection  to  the  surface 
of  the  particle  is  positive,  and  recalling  that  equilibrium  has  been 
neglected  in  the  gas  phase,  yields 


q 


II 

c 


AdT 

ar 


r 

P 


(21) 


where  the  second  part  of  Eq.  (21)  follows  from  the  assumption  of 
constant  gas  phase  specific  heat.  Taking  the  derivative  of  Eq.  (12) 
and  evaluating  the  resulting  expression  at  r  =  r  ,  yields 


^c  p  _  _ -exp(-K/Le) _ 

K.(h  -h  )  exp(-K/yLe)  -exp(-K/Le) 

CO  p 


(22) 


Given  K,  Eqs.  (19)  and  (20)  yield  the  four  element  mass  fractions 
at  the  surface  of  the  particle.  Equation  (11)  then  provides  four  of  the 
equations  to  determine  the  seven  unknown  species  mass  fractions, 
imposing  the  condition  of  chemical  equilibrium,  yields  the  three  more 
equations  needed  to  determine  the  Y..  The  following  dissociation 
reactions  were  employed 
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Since  the  1^  are  known  functions  of  temperature,  Eqs.  (4),  (19), 

(20)  and  (24)- (26),  allow  Y.  and  Y.  to  be  determined,  given  the 

XP  ZP 

ambient  composition,  K, and  tne  surface  temperature.  The  next  section 
will  consider  the  relationship  between  K  and  the  Y^  . 

P 

Constant  average  properties  were  employed  at  each  instant  of  time, 
using  methods  similar  to  earlier  work  in  this  laboratory  [9].  The 
methods  used  to  compute  properties  are  described  in  Appendix  A. 

3.2.2  Surface  Reactions 

The  reaction  of  carbon  is  a  complex  process,  particularly 
when  the  material  is  in  a  combustion  product  mixture  containing  many 
species.  Reaction  mechanisms  of  carbon  involving  0^ ,  CO^,  H^O,  H^,  OH 
and  0  can  play  a  role  [15,16].  During  the  present  analysis,  only  the 
first  three  species  0^,  CO^  and  H-O  were  considered.  The  present 
methodology,  however,  is  applicable  to  more  complex  treatment  of 
reaction  as  new  information  becomes  available. 


i 


16 


Other  problems  that  must  be  considered,  include:  1)  reaction  is 
not  just  limited  to  the  apparent  outer  surface,  since  pores  are  present; 
and  2)  the  presence  of  a  catalyst  influences  reaction  rates.  The  effect 
of  both  these  factors  was  treated  in  the  present  analysis  by  intro¬ 
ducing  empirical  area/reactivity  multiplication  factors  for  each 
reaction  contributing  to  the  consumption  of  carbon.  The  values  of  the 
area/reactivity  factors  were  specified  by  matching  predictions  and 
measurements. 

The  reaction  rates  for  carbon  with  oxygen  and  carbon  dioxide  em¬ 
ployed  the  same  kinetic  parameters  as  Libby  and  Blake  [13].  Tne 
reactions  are  assumed  to  proceed  by  the  following  global  expressions 


C  +  1/2  0, 


C  +  C02  2C0 


where  the  specific  rate  constants  are 


R1  “  K1  p02  ’  R2  ”  K2  PC02 


where 


K.  =  A.  exp(-E./RT)  ,  p.  =  pMY./M. 
iii  i  it 


The  specific  values  of  the  chemical  parameters  appearing  in  Eqs.  (29) 
and  (30)  are  summarized  in  Table  3. 

The  reaction  rate  expression  for  water  vapor  employs  the  results 
of  Johnstone,  et  al.  [17].  The  reaction  mechanism  is  as  follows  [13] 


17 


Table  3 


Summary  of  Reaction  Rate  Parameters 


i 

A 

E(kcal/gmol) 

i 

87100  kg/m^s  atm 

35.8 

2 

2470  kg/m^s  atm 

41.9 

3 

5.15  x 

10  ^  kg/ro^s  atm3 

32.7 

4 

9.42  x 

io-11 

-60.8 

5 

7.07  x 

IQ'16 

-79.3 

Slurry 

Reactivity  Factor,  a^*3 

noncatalyzed 

58.2 

catalyzed 

70.6 

2 

Assuming  a  surface  area  of  reaction  of  1.15  m  /g  of 
carbon,  which  is  the  average  surface  area  over 
the  period  of  reaction  [17]. 

^Reactivity  factor  taken  to  be  the  same  for 
reaction  with  0^,  CO^  and  f^O. 
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zz 


(31) 


C  +  H20  H2  +  C(0) 

C(0)  *  CO  (32) 


The  reaction  rate  of  carbon  due  to  the  presence  of  water  vapor  can  then 
be  expressed  as 


lSp«,o 

"* '  r^*V';'SV 

Z  Z 


(3  3) 


The  parameters  in  this  case,  appearing  in  Table  3,  were  obtained  from 
Ref.  17.  In  order  to  obtain  the  rate  constants  of  Table  3  for  these 
reactions,  it  was  necessary  to  estimate  the  actual  surface  area  for 
reaction  for  the  experiments  of  Kef.  17.  The  surface  areas  are  taken 
as  the  average  over  the  range  ot  reaction,  measured  by  nitrogen 
adsorpt ion  [ 1 7 | . 

The  total  reaction  ot  carbon  was  obtained  as  the  sum  of  the  rates 
or  02,  CO 2  and  HuO,  ignoring  potential  interactions  between  the 
reactants.  This  yields 


K 


JP 

I) 


3 

i  =  i 


(34) 


where  the  a j  appearing  in  Kqs.  (34)  are  t  lie  area/reactivity  multiplica¬ 
tion  factors.  These  parameters  were  varied  in  the  present  investigation 
in  order  to  match  burning  rate  results  to  predictions.  In  general  the 
presence  of  pores  and  a  catalyst  result  in  a  values  greater  than  unity 


19 


it 


r 


3.2.3  Agglomerate  Life-History  Computations 

Given  the  temperature  and  size  of  the  particle,  and  the 
^mbient  conditions,  solution  of  Eqs.  (4)— (31)  are  sufficient  to  provide 
Yi  ,  Yi  ,  q"  and  K.  The  variation  of  the  size  and  temperature  of  the 
P  P  c 

particle  are  determined  by  solving  equations  for  conservation  of 
particle  mass  and  energy. 

Conservation  of  mass  for  the  particle  is 


JE  <3  "rp3  p8>  * 


4iir 


(35) 


Assuming  that  the  particle  density  is  constant, 


dr 

-r-2-  =  ~  pDK/p  r 
dt  s  p 

introducing  K  from  Eq.  (16). 

Conservation  of  energy  at  the  surface  of  the  particle  yields 


(36) 


JE<T  V  ps  V  ■  4%2  <E*1  hi  +  "c  +  0?  (3,) 

where 

q"  =  O  e  (T  4  -  T  4)  (38) 

r  p  w 

and  T  is  the  temperature  of  the  enclosure  surrounding  the  hot  gas. 
Simplifying  Eq.  (34)  by  introducing  conservation  of  mass,  Eq .  (17), 
yields 


_E  = 

dt 


- — -~~r  [m"  (h 
(pscgrp)  p  cP 


hg)  +  q^  +  q';i 


I 

(39) 


During  some  experimental  conditions,  the  rate  of  radius  change 
of  the  particle  was  relatively  slow.  In  this  case,  the  particle  would 
reach  a  quasisteady  equilibrium  condition  where  its  rate  of  temperature 


change  with  time  was  small.  These  quasisteady  particle  temperatures 
were  determined  by  setting  the  left-hand  side  of  Eq.  (39)  equal  to 
zero  and  solving  the  equation  for  particle  temperature. 

Conservation  of  particle  momentum  would  yield  an  additional 
equation  to  determine  the  motion  of  the  particle.  Since  supported 
particles  were  considered  during  this  investigation,  these  results 
are  not  needed  and  will  not  be  considered  here. 

The  initial  conditions  for  Eqs.  (36)  and  (39)  are 


t  =  0, 


(40) 


Equations  (36)  and  (40)  were  numerically  integrated  on  the 
computer,  using  Gear's  method.  The  nonlinear  system  of  algebraic 
equations  which  must  be  satisfied  in  order  to  determine  V;  ,  ,  q" , 

P  P  c 

etc.,  were  solved  using  the  Newton-Raphson  method,  given  the  current 
values  of  T  and  r  ,  and  the  ambient  conditions,  at  each  instant  of 
time.  Quas?steady^burning  rates  were  determined  by  iterating  on  T  , 
until  the  left-hand  side  of  Eq.  (39)  was  zero.  ^ 


4.  Results  and  Discussion 


4.1  General  Observations 


Figure  4  is  an  illustration  of  the  variation  of  mean 
properties  along  the  centerline  of  the  propane  flame  [9].  Near  the 
injector  exit,  gas  velocities  and  the  concentration  of  propane  are 
relatively  high,  while  the  temperature  is  low  and  oxygen  is  absent. 
Moving  downstream,  the  gas  velocity  and  t he  concentration  of  propane 
decrease  monatonically .  The  temperature  increases  at  first,  reaching 
a  maximum  at  the  point  where  the  fuel  has  disappeared  and  com¬ 
bustion  product  concentrations  are  highest.  beyond  this  location 
of  the  flame,  the  flow  decays.  in  this  region,  temperatures  ami 
product  gas  concentrations  decrease,  while  oxygen  and  nitrogen  con¬ 
centrations  increase. 

Similar  to  our  previous  slurry  combustion  study  [4],  present  mea¬ 
surements  were  obtained  along  the  centerline  of  this  flame.  Also 
marked  on  Fig.  4  ate  various  combustion  regions  determined  for  both 
catalyzed  and  noncatalyzed  slurry  drops  during  our  earlier  work  14). 
Combustion  performance  for  each  of  these  regimes  described  in  fable  4. 
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Summary  of  Slurry  Drop  Combustion  Regions  in  the  Turbulent  Flame 


1 

I 


1 


The  present  measurements  were  confined  to  the  range  225  <  x/d  _  511). 

This  includes  the  full  combustion  (FC),  glowing  (G)  and  evaporation 
(E)  regions.  In  the  first  two  regions,  combustion  of  the  carbon  ag¬ 
glomerate  is  observed,  although  the  second  region  is  only  observed 
for  the  catalyzed  slurry.  In  the  evaporation  region,  combustion  of 
the  agglomerate  is  quenched.  This  region  is  of  marginal  interest  for 
the  present  study,  however,  it  is  desirable  to  predict  its  onset. 

Table  5  is  a  summary  of  flow  properties  at  each  of  the  drop  test 
locations. 

Figure  5  is  an  illustration  of  a  noncatalyzed  slurry  drop-life-  ’ 

iiistory  at  x/d  =  297.5,  where  agglomerate  reaction  rates  are  reasonably 
high.  Particle  diameter  and  temperature  are  plotted  as  a  function  of 

time.  The  time  scale  has  been  expanded  during  the  first  portion  of  i  ; 

the  plot,  when  the  liquid  is  gasified,  since  this  period  is  relatively 
short . 

The  four  main  stages  of  slurry  drop  combustion  are  quite  apparent 
in  Fig.  5.  The  initial  phase  of  the  process  involves  heat-up  of  the 
particle,  approaching  the  wet-bulb  temperature  for  liquid  evaporation 

(which  is  somewhat  below  the  boiling  temperature  of  the  liquid  [11]).  < 

Ln  this  period,  the  rate  of  liquid  evaporation  is  relatively  low  and  i 

the  particle  diameter  remains  nearly  constant.  As  the  surface  tem¬ 
perature  nears  the  wet-bulb  temperature,  the  rate  of  gasification 
increases,  beginning  the  second  stage  of  the  process.  During  the 
second  stage,  the  particle  diameter  decreases  while  the  particle 

temperature  remains  nearly  constant.  The  third  stage  involves  heat-up  ’ 

of  the  agglomerate,  once  all  the  liquid  has  evaporated.  The  variation 

of  particle  diameter  is  relatively  slow  in  this  period,  due  to  low 

carbon  reaction  rates  at  low  temperatures.  As  particle  temperatures 

increase,  reaction  rates  increase  as  well,  signaling  the  onset  of  the 

fourth  stage  of  the  process.  During  the  fourth  stage,  the  agglomerate 

reacts,  with  particle  temperatures  near  the  gas  temperature.  In  this 

regime,  the  particle  temperature  can  be  either  above  or  below  the  gas 

temperature,  depending  upon  the  magnitudes  of  energy  released  by  reaction 

and  heat  transport  by  convection  and  radiation.  These  quantities  tend  I 

to  vary  with  particle  size,  so  there  is  no  fixed  particle  temperature  j 

in  the  fourth  regime,  even  though  ambient  conditions  are  fixed.  As  j 

the  particle  becomes  small,  however,  convection  heat  transfer  dominates, 
and  the  particle  tends  to  approach  the  local  gas  temperature. 

It  is  evident  from  the  results  illustrated  in  Fig.  5,  that  heat-up  1 

and  reaction  of  the  agglomerate  comprises  the  bulk  of  the  lifetime  of 
the  particle. 

A . 2  Agglomerate  Burning  Kates 

Figure  6  is  an  illustration  of  predicted  and  measured  agglomerate 
burning  rates  for  noncatalyzed  drops.  Results  are  shown  for  various 
positions  in  the  range  297  £  x/d  <  510.  This  region  is  in  the  lean 
portion  of  the  propane  flame,  corresponding  to  an  equivalence  ratio  range 
of  0.761-0.272.  For  these  conditions,  reaction  of  carbon  with  oxygen 
dominates  the  process,  and  the  effects  of  CO.,  and  Ho0  are  small. 
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RATE  OF  CHANGE  OF  PARTICLE  DIAMETER  (pm/s) 


In  Fig.  b,  predictions  have  been  matched  to  the  measurements  at 
x/d  =  382.5,  by  selecting  a  reactivity  factor  of  a.  =  58.2  for  all 
the  reactions.  This  parameter  was  then  held  constant  for  the  re¬ 
maining  computations,  a  procedure  which  yields  a  reasonably  good 
correlation  of  the  remaining  measurements,  including  both 
effects  of  particle  size  (in  the  range  400-730  pm)  and  position  in 
the  flame.  The  low  reaction  rates  indicative  of  the  evaporation 
region  were  also  predicted  reasonably  well.  Thus  the  model  demon¬ 
strates  some  potential  for  predicting  quenching  of  agglomerate 
reaction. 

There  are  two  main  combustion  limits  for  carbon  agglomerates. 

The  first  limit  involves  reaction  rate  control,  where  reaction  rates 
are  relatively  low  in  comparison  to  diffusion  rates, and  the  concen¬ 
trations  of  gaseous  reactants  at  the  surface  of  the  particle  are 
nearly  the  same  as  in  the  surroundings.  This  regime  is  characterized 
by  the  burning  rate  being  relatively  independent  of  particle  diameter, 
e.g.,  results  for  x/d  382.5  in  Fig.  6  approach  the  reaction  rate 
limit . 

The  second  limit  involves  diffusion  control,  where  reaction  rates 
are  high  in  comparison  to  diffusion  rates,  and  the  concentrations  of 
reactant  gases  at  the  surface  of  the  particle  are  small.  The  variation 
of  the  burning  rate  is  influenced  by  particle  diameter  in  this  region, 
showing  a  dependence 


n,  i/r  <0-5  to  X-0> 

P 


(41) 


The  power  of  r  ,  in  Eq.  (41),  tends  to  increase  as  the  Reynolds  number 
of  the  particlE  decreases.  From  this  argument,  it  is  apparent  that 
the  higher  reaction  rates  illustrated  in  Fig.  6  (x/d  =  297.5,340.0)  are 
approaching  the  diffusion  controlled  limit.  Diffusion  rates  become 
large  for  small  particles;  therefore,  all  regions  of  the  flow  would 
eventually  exhibit  kinetic  control  for  sufficiently  small  particles. 

Figure  7  is  an  illustration  of  predicted  and  measured  burning  rates 
for  a  catalyzed  agglomerate.  These  test  results  were  limited  to 
x/d  =  297.5,  340.0  and  382.5.  Reaction  of  the  agglomerate  is  nearly 
diffusion  controlled  at  the  first  two  positions.  Therefore,  burning 
rates  with  and  without  catalyst  are  essentially  the  same  at  these 
locations.  At  x/d  =  382.5,  however,  the  burning  rate  is  more  nearly 
reaction  rate  controlled.  At  this  position,  the  use  of  the  catalyst 
results  in  a  20-25  percent  increase  in  the  burning  rate.  Naturally, 
the  effect  of  the  catalyst  is  even  larger  at  positions  farther  down¬ 
stream,  e.g.,  adding  catalyst  reduces  the  lean  limit  for  quenching  from 
0.37  to  0.21  [4] . 
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The  presence  of  catalyst  is  treated  in  the  theory  by  adjusting  the 
area/reactivity  factor.  Selecting  a.  =  70.6,  for  all  the  reactions, 
yields  the  theoretical  results  illustrated  in  Fig.  7.  With  this  choice, 
the  predictions  are  in  reasonably  good  agreement  with  the  measurements. 

Another  property  of  catalyzed  slurries  was  their  tendency  to 
spontaneously  shatter  from  time  to  time.  This  shattering  occurred 
after  all  the  liquid  had  apparently  evaporated.  This  behavior  could 
be  influenced  by  the  presence  of  the  thermocouple  support,  and  was  not 
consistently  observed.  If  shattering  does  persist  for  small  unsup¬ 
ported  particles,  this  would  provide  another  (physical)  mechanism  for 
the  improved  combustion  properties  observed  for  catalyzed  slurries 
[3].  Tests  with  unsupported  particles  having  sizes  more  typical  of 
sprays  will  be  required  to  resolve  this  question. 

4 . 3  Agglomerate-Life-Histories 

Agglomerate-life-histories  were  computed  at  x/d  =  255-510, 
for  noncatalyzed  slurries.  The  resulting  comparison  between  predicted 
and  measured  particle  diameters  and  temperatures  as  a  function  of  time 
are  illustrated  in  Figs.  8-13.  For  reference  purposes,  the  local  tem¬ 
perature  of  the  flow  is  also  shown  on  each  plot. 

The  initial  condition  for  the  computations  illustrated  in  Figs. 

8-13,  was  the  temperature  and  diameter  of  the  agglomerate,  when  all 
the  liquid  had  evaporated.  The  comparison  between  predictions  and 
measurements  : s  seen  to  be  quite  good.  This  is  partly  expected,  since 
the  model  was  capable  of  correlating  the  burning  rate  measurements 
for  this  range  of  flame  conditions.  However,  the  results  show  that 
the  model  gives  a  good  representation  of  the  heat-up  of  the  particle 
as  well,  where  reaction  effects  are  small,  for  various  particle  sizes 
and  ambient  temperatures,  compositions  and  velocities,  c.f.  Table  5. 

At  higher  temperature  levels,  heat-up  is  also  influenced  by  reaction 
and  radiation  to  the  surroundings.  It  is  promising  that  the  model 
appears  to  be  capable  of  treating  these  effects  as  well. 

The  range  of  conditions  illustrated  in  Figs.  8-13  is  such  that 
particle  heat-up  is  terminated  at  temperatures  both  above  and  below 
the  local  gas  temperature.  This  behavior  is  due  to  the  opposing 
effects  of  reaction,  tending  to  increase  particle  temperatures  above 
the  gas  temperature,  and  radiation,  which  tends  to  reduce  the  particle 
temperatures  below  the  gas  temperature.  In  all  cases,  however,  con¬ 
vection  would  dominate  the  process  as  the  particle  becomes  small  and 
the  particle  would  finally  tend  to  approach  the  gas  temperature.  This 
is  not  observed  in  Figs.  8-13  due  to  the  finite-sized  thermocouple  bead 
used  to  mount  the  slurry  drop. 

4 . 4  Limitations  of  the  Study 

While  the  model  shows  an  encouraging  capability  to  predict  the 
combustion  properties  of  carbon  agglomerates,  several  limitations  of  the 
present  evaluation  should  be  recognized.  The  main  limitation 
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TIME  AFTER  ALL  LIQUID  EVAPORATED  (s) 
Figure  8.  Agglomerate-life-history  at  x/d  =255.0 
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TIME  AFTER  ALL  LIQUID  EVAPORATED  (s) 

Figure  10.  Agglomerate-1 i f e-history  at  x/d  =  34 
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10.  20.  30.  40.  bU. 

TIME  AFTER  ALL  LIQUID  EVAPORATED  (s) 

Figure  12.  Agglomerale-1 i fe-history  at  x/d  =  4f >7.3. 
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EMPERATURE  (k) 


of  this  study  is  the  large  particle  sizes  examined, 

400-1000  ym,  in  comparison  to  the  sizes  expected  in 
practical  combustors.  The  major  difficulty  involves  the  present 
specification  of  area/reactivity  multiplication  factors — which  were 
taken  to  be  constant  during  the  present  study.  For  large  particles, 
diffusion  of  reactant  gases  is  not  expected  to  penetrate  very  far 
into  the  particle,  while  very  small  particles  would  have  reactant 
gases  penetrated  throughout  their  volume.  Therefore,  the  relative 
proportions  of  apparent  surface  to  actual  reaction  areas  will  tend 
to  increase  as  particle  sizes  become  small.  Therefore,  it  is  ex¬ 
pected  that  the  area/reactivity  multiplication  factor  would  tend  to 
increase  as  particle  sizes  are  reduced.  The  magnitude  of  this  change 
can  only  be  established  by  additional  study.  Naturally,  this  effect 
would  only  be  important  outside  the  diffusion  controlled  limit. 

Present  test  conditions  were  not  sufficiently  varied  to  establish 
the  appropriate  area/reactivity  multiplication  factors  for  the  carbon 
dioxide  and  water  vapor  reactions,  since  oxygen  dominated  carbon 
reaction  when  the  basic  reaction  parameters  of  Table  5  are  used  in  the 
analysis.  It  is  unlikely  that  the  effect  of  pores,  surface  properties 
or  catalyst  would  act  the  same  for  each  reactant.  Therefore,  use  of 
the  present  findings  should  be  limited  to  conditions  where  reaction  of 
carbon  with  oxygen  dominates  the  process.  Possible  interference  effects 
between  the  various  reactants,  and  whether  all  the  appropriate  reactants 
were  considered  also  remains  to  be  resolved. 

Another  variable  of  significant  practical  importance  is  the  effect 
of  pressure.  The  present  results  were  limited  to  atmospheric  pressure 
and  whether  the  reaction  expressions  used  here  yield  the  proper  pressure 
dependence,  must  still  be  evaluated. 


}7 


-  * 


Appendix  A 


Computation  of  Properties 
A. 1  Mean  Property  State 

The  average  properties  of  the  gas  surrounding  the  droplet  were 
computed  at  the  mean  state,  defined  as  follows: 


0  =  a0  +  (1-a)  0 

avg  p  °° 


(A-l) 


where  0  is  a  generic  property  at  a  specific  temperature  and  mass 
fraction,  a  can  range  between  0  and  1  and  a  value  of  0.9  was  selected 
for  the  present  investigation  as  suggested  by  Faeth  and  coworkers 
[9,10]. 


A.  2  Mixture  Thermal  Conductivity 

The  thermal  conductivity  of  the  gaseous  mixture  was  calculated 
with  the  Mason  and  Saxena  formulation  of  the  Wassiljewa  equation  [18]. 
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x .  A . 
1  1 


f  X.4>.. 

j-1  J  1J 


(A- 2) 


where 


ij 


aU2  U  +  Mj> 


-1/2 


P.  1/2  M.  1/4 

1  *  <7T>  ‘S1’ 

J  1  -J 


(A -i) 


The  thermal  conductivity  and  viscosity  for  gases  other  than  JP-10 
was  obtained  from  Ref.  19.  The  thermal  conductivity  of  JI-10  was  de¬ 
termined  using  the  Miscic  and  Thodos  method  [18]  and  the  viscosity  of 
JP-10  was  determined  using  the  Reichenberg  method  118].  Equation  (A-2) 
was  evaluated  at  the  average  temperature  and  composition  of  the  gas 
phase.  In  Eq.  (A-3),  ]>„  is  set  equal  to  1,  when  i  -  j  . 


Mixture  Viscosity 


A. 3 


The  viscosity  calculations  for  the  mixture  employed  the  method 
of  Wilke,  cited  in  Ref.  18. 

For  a  mixture  of  N  components 

VA-4 ) 

The  viscosity  data  used  in  the  calcu- 
manner  as  described  in  Section  A. 1. 
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x .  f  .  . 
J  ij 


where  <$> .  .  is  given  by  Eq.  (A-3). 
lations^as  obtained  in  the  same 


A. 4  Mixture  Diffusivity 

When  all  of  the  liquid  fuel  has  evaporated  and  only  a  carbon 
agglomerated  particle  remains  the  diffusion  coefficient  of  the  gas  was 
taken  for  simplicity  as  the  binary  diffusion  coefficient  of  oxygen 
through  nitrogen 


D 


°2’N2 


(A-r>) 


Use  of  other  species  combinations  yielded  di f f usivit ies  of  similar 
magnitude . 
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